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SUMMARY 


A one-dimensional computer model of the Earth’s troposphere and strato- 
sphere is described. The model includes a photochemical reaction set of 
87 kinetic processes and 28 photolytic reactions^ with vertical transport 
being taken into account by an eddy diffusion parameter. Temperature pro- 
files are calculated with a radiative-convective code which balances infra- 
red atmospheric emission with solar absorption. The temperatures computed 
may be included in the iterative solution procedure; thus, temperature- 
chemistry coupling effects may be investigated. An anisotropic multiple- 
scattering code may also be included as a part of the coupled calculation. 

The model may be used in either a direct insolation or diurnal average 
computation. The diurnal averaging procedures are described fully. A com- 
plete listing of the photochemical mechanism is given, along with current 
values for rates and cross sections. Boundary conditions for a sample cal- 
culation and the resulting vertical profiles of species concentrations and 
temperatures are given. It is concluded that the model gives a reasonable 
representation of many physical and chemical aspects of the globally averaged 
atmosphere. 


INTRODUCTION 

Study of the atmosphere is inherently complicated, involving dynamical, 
chemical, and radiative processes as well as their interactions. To develop 
a mathematical description of the atmosphere, it is necessary to understand 
individual processes. For this reason, atmospheric scientists have developed 
a hierarchy of models, each serving a specific purpose and typically concen- 
trating on one aspect of the overall problem while treating other aspects 
in a simplified manner. One-dimensional models, which represent globally 
averaged conditions, were designed primarily to study chemical processes. 
Species transport is included by using an eddy diffusion parameterization 
that gives the net effect of several dynamical processes (mean flow, plan- 
etary waves, etc.). Three-dimensional models, on the other hand, provide 
detailed treatment of dynamics, while either ignoring chemistry completely 
or employing very truncated chemical schemes. Despite their relatively 
simplistic treatment of transport, one-dimensional models have performed 
remarkably well in predicting the mean vertical distributions of strato- 
spheric trace species when compared with observed profiles. Because of this 
fact, one-dimensional models have been used extensively in the past (see, 
for example, refs. 1 through 7) and, because of their computational effi- 
ciency, they will continue to be used widely to examine new kinetic mech- 
anisms before incorporating them into more complicated models and to assess 
the impact of man-induced perturbations on the physical state of the 
atmosphere. 



A one-dimensional steady-state radiative-convective-photochemical model of 
the atmosphere has been developed recently at the Langley Research Center to 
study the chemical and thermal structure of the troposphere and stratosphere. 

The complete chemical model extends from ground level to the approximate height 
of the stratopause (0-55 km) . Although many questions of current interest are 
nonsteady (for example, the release of chlorofluorome thanes) , the model was not 
developed to predict the state of the stratosphere, but rather was developed to 
investigate atmospheric phenomena, specifically coupling between thermal and 
chemical processes. For this purpose, a steady-state model was found to be 
sufficient. In addition to chemical- thermal feedback, the model also incorpo- 
rates molecular scattering processes and reflection of radiation from the 
ground and clouds. Studies may be made of the effects of injections of nitrogen 
oxides, water vapor, carbon dioxide, carbon monoxide, and chlorofluorome thanes, 
as well as investigations of the impact of changes in solar radiation on the 
atmospheric thermal structure. The model is fully coupled in the sense that 
the calculated thermal structure is consistent with the vertical profiles of 
the chemical species. 

This paper presents the Langley model of the atmosphere. The following 
features of the model are relatively completely, but concisely, described: 

System of chemical and photolytic processes 
Atmospheric radiative transfer processes 
Computation and solution methods 
Current input data 

Representative examples of calculated results 


SYMBOLS 

infrared band absorptance, m“^ 

specific heat for constant pressure, J-kg“^-K“^ 

Planck blackbody function, W-m“^ 

see equation (17), W-m"^ 

radiative flux, W-m”^ 

fractional cloud cover 

average pressure scale height, km 

mean intensity of radiation, W-m"^ 

photolysis rate for ^>th photolytic process, s“^ 

Jacobian matrix (see eqs. (16)), s“^ 

vertical eddy diffusion coefficient, m^-s"^ 


A 

ex 

f(K,K) 

Fr 

^c 

H 

I 

J 

Kz 
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II 


k 

kjl 

L 

n 

P 

P 

Q 

q 

r 

S 

T 

t 

U 

V 


Boltzmann constant, 

chemical rate coefficient for Jlth kinetic process, m^-s^^-particle”*^ 
or m^-s"^ -particle'^ 

chemical loss rate, s“^ 

number density, particles-m”^ 

chemical production rate, particles-m“^-s~^ 

pressure. Pa 

heating rate, K-s“^ 

precipitation lifetime, s 

ratio of nighttime to daytime average concentrations of chemical 
species 

net chemical production of chemical species, particles-m“^-s“^ 
temperature, K 
time, s 

length of day, s 
length of night, s 

length of total diurnal cycle, 24 hr 

unit matrix 

velocity vector , m-s“^ 


z altitude, km 

a solar zenith angle, rad 

diffusion-related quantities (see eq. (13)), s”^ 

Y see equation (11), km~^ 

6 defined in equation (28) 

ri quantum efficiency for photolysis 

X wavelength of radiation, m”^ 

^ dimensionless optical path length 
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P gas density, kg-m*"^ 

0 absorption cross section, 

T dimensionless optical depth 

$ particle flux, particles~m”2-s"^ 

X mole fraction 

Subscr ipts: 

a light-absorption process 

c value with inclusion of clouds 

1 chemical species 

j nodal point in numerical calculation 

a chemical or photolytic process 

nc clear sky condition 

o surface value 

s light-scattering process 

T value at tropopause 

°° value at top of atmosphere 

Super scr ipts : 

d daytime value 

n nighttime value 

°° total atmospheric optical depth 

An arrow over a symbol denotes a vector. Brackets around a chemical species 
(e.g. fCl] ) denote its number density. 

GOVERNING EQUATIONS 

A pictorial representation of the physical and chemical processes con- 
sidered is given in figure 1 . The species concentrations are governed by the 
standard conservation equation, 


4 


9t 


+ V • iniV 


dt 


chem 




= P-i 


- 


(1 ) 


which in the one -dimensional approximation becomes 


dnj_ 9$i 

+ = ^i - (2) 

at 3z 


Here the chemical source term is expressed as the difference between produc- 
tion Pi and loss Lini of the ith species. The quantity is the verti- 

cal flux which in a one -dimensional model is parameterized as being propor- 
tional to the concentration gradient 


$i 


-nKz 


azln ^ 


(3) 


where n is the total number density and is the vertical diffusion coef- 

ficient and represents the effect of large-scale processes on the vertical 
transport of species, ^en the equation of state (n = p/kT) and an average 
pressure scale height H are introduced, equation (3) may be rewritten as 


^i 





(4) 


which is one of the more common forms of the flux expression used in one- 
dimensional models. 


In the steady-state solution of the conservation equation, the time deriv- 
ative of n^ vanishes; the equations to be solved for the species concentra- 
tions are then 



(5) 


The equations are generally coupled because the chemical production and loss 
terms depend on other species as well as on temperature and altitude. The 
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explicit finite-difference equations used in solving the species conservation 
equations are described in detail later. 

The production and loss terms in equation (5) are calculated from speci- 
fied reaction rates. For photodissociation processes, such as AB + hv A + B, 
the rates depend on the local solar intensity in addition to the species 

number densities and spectral absorption coefficients a(X) : 


Jjl(z) 


■I 


dX a5,(X) ri£(X) 




-r I - 


(X) ni sec a dz' 


(6) 


where n j, is the quantum efficiency, is the mean solar intensity at the 

top of the model atmosphere, and a is the solar zenith angle. 

Equation (6) may be used to compute the photolysis rates for a fixed sun 
(a = Constant) or over a diurnal cycle by using the appropriate time dependence 
of a. The diurnal calculation incorporates several approximations for averag- 
ing of the nonlinear product terms in the species conservation equation; these 
approximations are described more completely in a subsequent section. 

The photodissociation rates are also affected by light-scattering pro- 
cesses in the atmosphere. When these effects are included, scattering factors 
Js,5,/Ja,il calculated from the following relation: 






Ico(^) exp/[xs(X) - Ts(X/Z) + Ta(X) 



I 


dX + J I(X,z) dX 


(7) 


where I is the diffuse intensity integrated over all solid angles, Xg 

oo 

and Xg are the total optical depths due to scattering and absorption, and Xs 
and Xa are the optical depths due to scattering and absorption at the altitude 
of interest. The quantity given by equation (6), is the photolysis rate 

in an atmosphere which only absorbs. The first term of equation (7) represents 
the contribution to photolysis of the direct beam attenuated by scattering and 
absorption, while the second term is the contribution of radiation which has 
been scattered one or more times. 

If the vertical temperature profile is not specified, the energy conserva- 
tion equation must be introduced. In the Langley model, the vertical tempera- 


6 



ture structure is calculated with a radiative-convective code wherein only radi- 
ative transfer processes are taken into account. In this situation, the energy 
equation reduces to 


Q = 


3T 

3t 


1 

V • 

pcp 


Fr 


(8) 


where Fj^ is the radiative flux vector which includes the absorption of solar 
radiation and the absorption and emission of infrared radiation by the optically 
active constituents, O 3 , H 2 O, and C 02 * The absorption bands of these species 
are included by using a band absorptance and emissivity formulation, putlined 
in reference 8 , which avoids the time-consuming and tedious integrations with 
respect to frequency over several thousands of spectral lines. 

Because the stratosphere is nearly in a state of radiative equilibrium, 
temperatures in this region are calculated by setting Q equal to 0 in equa- 

-y 

tion ( 8 ) . The temperature-dependent terms in V * Fr are used to determine 
the mean thermal structure. The solar heating term in the flux divergence is 
multiplied by 1/2 to account for the absence of heating during the night. With- 
in the troposphere, the temperature profile is taken to decrease with height in 
accordance with the mean observed tropospheric lapse rate, which approximates 
energy transport by convective processes and latent energy release by water 
vapor within the troposphere. Further details concerning the temperature cal- 
culation and coupling the radiative and photochemical codes are given in the 
next section. 


SOLUTION METHODOLOGY 

In this section, the procedures used to obtain solutions to the equations 
given in the preceding sections are described. To facilitate the description 
of the various aspects of the calculation, this section is divided into the 
following subsections: 

Solution to the Photochemical-Transport Equations 
Temperature Calculations 

Coupled Calculations With Temperature and Scattering- 
Albedo Effects 
Diurnal Calculations 


Solution to the Photochemical-Transport Equations 

As mentioned in an earlier section, the Langley model has been used to 
obtain steady-state solutions for the species distributions. For this case, 
the conservation equation for species i reduces to 
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+ Si = 0 


(9) 




dz 


where Si = Pi - niLi is the net photochemical production or loss of the 
ith species. Equation (9) is solved between the ground and 55 km by divid- 
ing this region into N subintervals. For the present calculations, N = 40 
has been used. Special care is taken to insure that a computational node is 
placed at those points where the vertical diffusion coefficient undergoes a 
step discontinuity. 

A flux-conserving finite-difference representation of equation (9) is 
obtained by integrating at the jth node between Zj - Azj_-|/2 and zj + Azj/2 
to give 


J -»Zj+Azj/2 

Si(z') dz* = 0 (10) 

z j— Az j_i /2 

where A 2 j = Zj+i - Zj and the subscript j + 1/2 refers to the midpoint 
between the (j + 1 ) th and jth nodes. The flux (eq. (4)) is evaluated at the 
midpoint between two adjacent nodes with a second-order, central-difference 
approximation for the first derivative, which yields 


^i,j+l/2 


“*^z,j+l/2 


ni,j+i - ni,j 


Azj 


1 

+ - Yj+i/2(ni,j+i 


+ 



( 11 ) 


1 dT 1 

where y = + — . An expression analogous to equation (11) can also be 

T dz H 

written for ^i,j-i/2* The integral term in equation (10) is approximated as 


I 


Zj+Az j/2 


Zj-Azj_i/2 


Si(z') dz' 


^(Azj_i+ Azj)Si,j 


( 12 ) 


Substitution of equations (12), (11), and the corresponding expression for 
into equation (10) and rearrangement of terms yield 

“^l/j"i/j-l * ^2,j*^i,j ■ ^3,j«i,j+i - Si^j = 0 (13) 

where the 3 's are diffusion-related quantities involving the spatial step 
size. In the Langley model, each species distribution is determined from its 


8 



own conservation equation, rather than the species being grouped into families 

and the families tracked. One can therefore define a vector nj whose compo- 
nents represent all the unknown species concentrations at the nodal point j, 
and then equation (13) can be rewritten in matrix form as 


-3 





+ 3 


2, j 




j+1 



( 14 ) 


where U represents the unity matrix. Equations (1 3) and (1 4) are valid only 
for the interior nodes, but at the upper and lower boundaries similar expres- 
sions can be developed. For example, if number density is specified at the 
lower boundary, then the first term on the left side of equation (13) is a known 
quantity which can be transposed to the right side to give 


- Sij = Si^ini^O 


( 15 ) 


where j = 0 denotes the lower boundary. When flux is given, an expression of 
the same form as equation (15) results except that the right side now contains 
the specified flux. This latter expression and equation (1 5) can also be put 
into matrix form. 

Equations (1 4) represent a set of equations which are nonlinear and coupled 
because of chemical interactions among the various species. Further nonlinear- 
ities arise if interactions between the temperature field and species distribu- 
tions are allowed. For a prescribed temperature profile, equations (14) are 
solved by a Newton-Raphson iteration. Linearization of equations (1 4) about 

an initial species distribution nj yields the following for the correction 

terms 

Bl,jU 6nj_i + (B2,jU - J) 6nj - Ba^jU 6nj+i = ~Rj(nj) 6) 

where Rj represents the residual of equations (1 4) . The quantity J repre- 
sents the Jacobian matrix of the net chemical term at node j; that is. 


J 


if 


3Si 



Equations (16), along with the corresponding equations for the lower and upper 
boundaries, form a tridiagonal block system which can be solved efficiently 


for 6nj (see, for example, ref. 9) . The iterations are continued until the 
convergence criterion is satisfied. For the Langley model, the convergence 
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criterion is |e| ^ 0 . 001 , where e is the maximum relative change in number 
density between two successive iterations over all computational nodes and all 
species to which equation (9) is applied. The number of iterations needed for 
a solution depends on how close the initial distribution is to the final result, 
but typically 15 to 25 iterations are required. 


Temperature Calculations 

Given distributions of the radiatively active constituents - ozone, water 
vapor, and carbon dioxide - and an appropriate cloud model, one can compute a 
globally averaged vertical temperature structure by using a radiative-convective 
model. Since the radiative-convective code employed in the Langley model is the 
one developed in reference 8 and is described fully there, only a brief descrip- 
tion of its essential features is presented here. For further details, the 
reader is referred to the original paper. 


Stratospheric temperatures are calculated by assuming that the stratosphere 
is in radiative equilibrium; that is, the net radiative heating Q is zero. 

The term Q includes both longwave and solar contributions due to the various 
absorption bands of O 3 , H 2 O, and CO 2 . Within the troposphere, the temperature 
is assumed to decrease with height in accordance with the observed globally 
averaged lapse rate of -6.5 K/km. The temperature calculations then reduce to 
finding the surface temperature, the tropopause height, and the stratospheric 
temperature profile. Surface temperature is determined from the condition 
that, at the top of the atmosphere, the net outgoing longwave flux balances the 
net solar flux; that is, the net radiative flux Pr = 0. 


The stratospheric temperature distribution is given by the solution to the 
nonlinear, integral expression for the net radiative flux divergence. With 
the band absorptance formulation of reference 8 , the flux divergence for a 
particular absorber can be written generally as 


<3Fr(C) 

~dT" 


dA(C) nK 

-ex( 0 ) — - \ F(C,C) 

d^ 


rKi pCo 

J F(C,C’) dC* -J F(?,^’) d^' (17) 


where 


F (?,0 


dex(K’) (dA)(lC - 5'|) 
dC ^ 
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The quantity ^ is the dimensionless optical path length, which is measured 
fran the top of the atmosphere downward, and A is the band absorptance. 
Particular forms for ^ and A for the various bands of O 3 , H 2 O, and CO 2 
are given in reference 8 . The quantities and refer to the value 

of ^ at the tropopause and at the surface, respectively. The parameter 
is the Planck function evaluated at the band center. 

Since equation (17) is an integral equation, its solution is computation- 
ally time consuming and tedious. However, in reference 8 , it is shown that 
equation (17) may be closely approximated by the following expression: 


< 3 ? 


-ex (?) 


dA(?) 

< 3 ? 



F(?,?*) d^* 



F(?,?*) d^» 


(18) 


The only term in equation (17) that is neglected in equation (18) is the 
exchange of energy between the level ? and the region of the stratosphere 
lying above this level. This exchange of energy between the level ? and the 
upper region makes a small contribution when compared with the terms retained, 
since the absorber amount and the atmospheric pressure decrease exponentially 
with altitude. 

Within the stratosphere, the temperature distribution is obtained by let- 
ting Q = 0 (i.e., radiative equilibrium) where 


Q = QH20 + QCO2 Q03 Qs 


(19) 


The quantities QH 2 O' ^C 02 ' ^03 denote the longwave contributions to the 

radiative heating by water vapor, carbon dioxide, and ozone and are given by 
equations analogous to equation (18); Qg denotes the solar heating due to 
absorption by H 2 O, CO 2 , and O 3 . Use of equation (18) greatly simplifies find- 
ing the temperature distribution since equation (19) can new be treated as a 
nonlinear, algebraic equation, provided that the temperature at the lowest 
level within the stratosphere is obtained first, and then the temperature at 
the next higher level, and so on. The integral terms in equation (18) repre- 
sent contributions frm the atmospheric regions belcw the level ?, which are 
easily evaluated if the previous procedure is followed. At each computational 
level, the temperature is obtained from equation (19) by employing a Newton- 
Raphson iteration technique. The temperature is considered to have converged 
when the temperature change between two successive iterations is less than 
0.1 K. 


As mentioned earlier, the tropopause height must also be determined. This 
height is taken to be the highest level within the troposphere for which the 
lapse rate in the layer between this level and the next higher level becomes 
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dT 

smaller than the specified tropospheric lapse rate; that is^ — + 6*5 « 0. 

dz 

This level can be found by advancing upward from an initially assumed tropo- 
pause level which lies below the expected tropopause height. Temperatures are 
calculated by solving equation (1 9) . As long as the lapse rate within the layer 
above the assumed tropopause level exceeds the specified tropospheric lapse rate, 
the lapse rate is set equal to the prescribed value; this procedure is repeated, 

dT 

advancing from level to level, until the condition — + 6.5 « 0 is satisfied, 

dz 

Thus, for a specified surface temperature, the temperature distribution is cal- 
culated by first computing the tropopause height and then calculating strato- 
spheric temperatures beginning at the level just above the computed tropopause, 

Calculation of the surface temperature involves the following steps: 

( 1 ) assume a surface temperature; ( 2 ) compute the tropopause height and strato- 
spheric temperature profile by the method just explained; (3) calculate the net 
radiative flux at the top of the atmosphere; and (4) calculate a new surface 
temperature by a Newton-Raphson technique given by 


To (new) = To(old) 


Fr 

dFR/dTo 


( 20 ) 


Steps 2 to 4 are repeated until the condition = 0, which represents globally 
averaged conditions, is satisfied. Further elaboration on the procedures used 
for calculation of the temperature structure is given in reference 8 . 


Coupled Calculations With Temperature and Scattering-Albedo Effects 

Figure 2 illustrates schematically the procedure employed to include inter- 
actions between the temperature field and species concentrations. As can be 
seen, the photochemical and radiative codes are maintained as separate entities, 
with required information (namely, species profiles, primarily O 3 and H 2 O, and 
temperature distributions) being passed between the two program units. In addi- 
tion, the temperature calculations are lagged one step behind the calculations 
of species concentrations. To ensure that the concentrations and temperatures 
are consistent with one another, the cycle (fig. 2 ) is continued until the tem- 
peratures from two successive profiles differ by less than 0.25 K at all alti- 
tude levels. As long as temperature coupling is moderate, this procedure works 
satisfactorily. For stronger coupling, the number of iterations and the com- 
putation time increase. Typically, five to seven iterations on the temperature 
distribution are needed for a fully converged solution. 

Molecular scattering from air molecules and solar reflection from the 
surface and clouds are included in the calculation of species profiles accord- 
ing to the following procedure. For each photolysis process, the altitude 
variation of the ratio Js/'^a (see eq. (7)) is computed initially with a 
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multiple-scattering code described in this section, where Jg denotes the 
photolysis rate when molecular scattering and albedo effects are included and 
Jg is the photolysis rate without these effects, that* is, the photodissocia- 
tion rate in a purely absorbing atmosphere. To incorporate scattering and 
albedo effects, this array of ratios is stored and multiplied by the values of 
Jg which are calculated in the photochemical code. Because of the computa- 
tional demands of the multiple-scattering code, the ratio Js/»^a held fixed 
and is not reevaluated, as is done with the temperatures. This procedure is 
based on the assumption that Jg/Ja remains invariant to changes in the species 
distributions, particularly those which absorb at solar wavelengths. 

The value of Jg/'^a represents the weighted mean between two calculations, 
one for clear sky conditions and one with clouds included; that is. 


Js/Ja = n - fcM>^s/Ja)nc + ^c(^s/^a)c (21) 


where (Js/*^a) c denotes the ratio for full cloud cover, (d's/'^a^nc 
value for clear sky conditions, and f^;, represents the fractional cloud cover, 
which is taken as 0.5. The value of (Js/'^a^c computed for a single cloud 
layer at some effective altitude, typically taken as 5 km. 

The absorption and scattering properties of the atmosphere were calculated 
from the following data: total number densities of O 2 and N 2 from reference 10 , 

O 3 number density from the mid-latitude summer distribution of reference 11 (or 
alternatively from a photochemical calculation without scattering) , spectral 
intervals and cross sections of O 2 and O 3 from reference 12, solar flux 
from references 12 and 13, and the Rayleigh scattering cross sections from 
reference 14. 

An anisotropic multiple scattering radiation transfer code (ref. 15) has 
been employed to calculate the diffuse intensity field, I(X,z), at 40 equally 
spaced intervals between 0 and 50 km and for 110 frequency intervals between 
0.175 ^m and 0.735 urn. This field used in conjunction with equation (7) permits 
the evaluation of the scattering factors for the photodissociation processes of 
interest. 


Diurnal Calculations 

Species whose chemical lifetimes are short, less than a few hours, undergo 
significant variations near sunrise and sunset. Even species with lifetimes 
of 1/2 day, such as N 2 O 5 , have appreciable daily variations (ref. 16). Such 
variations can be important in determining the impact of perturbations on the 
ozone concentration for two reasons. First, use of a constant, average sun can 
result in underestimation of the true average concentration of a species which 
is particularly serious if that species is a sink for more reactive species. 

An example of this is N 2 O 5 . Calculations in reference 16 indicate that use of 
a 24-hour average sun underpredicts the average concentration of N 2 O 5 near 
25 km by roughly 2 orders of magnitude compared with results which take into 
account its diurnal variation. For a given total odd nitrogen concentration. 
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more N 2 O 5 implies smaller concentrations of NO and NO 2 to destroy ozone. A 
second reason to account for diurnal variations in species concentrations is 
accurate calculation of nonlinear terms occurring in the loss expression for 
odd oxygen. For example, the rate of destruction of odd oxygen by chlorine is 
given by 


Lci = 2k45[C10] [0] 


( 22 ) 


where k 45 is the rate constant for the reaction 


(45) CIO + 0 Cl + O 2 


The concentration of atomic oxygen decreases rapidly at sunset by reaction with 
O 2 to form O 3 , and when CIONO 2 is included in the chlorine chemical chain, the 
CIO concentration also decreases rapidly at sunset for altitudes below about 
40 km because of reaction with NO 2 to produce CIONO 2 at night. Thus, the aver- 
age destruction rate over 24 hours can be written as 


^C1 


2k45 


2 


(23) 


where superscript d represents the daytime average- However, when a 24-hour 
average sun is used, the destruction rate is evaluated as 


Lci = 2k45 



(24) 


with the bar denoting the 24-hour average, which, for CIO and 0, is related to 

their daytime average value according to = - n^^. Substitution of this last 

result into equation (24) and comparison with equation (23) show that use of a 
24-hour average sun causes the destruction rate of odd oxygen due to chlorine to 
be underestimated by approximately a factor of 2. This simple estimate agrees 
reasonably well with the more exact calculations of reference 17, indicating 
that the factor is 1.6 at 30 km. 

To accurately account for the impact of man- induced perturbations on ozone, 
the diurnal effects discussed in the previous two paragraphs must be included in 
the calculations. However, a direct time integration of the fully time-dependent 
equations is impractical because of excessive demands on computer resources when 
model simulation times of several decades are required, as in the chlorine-ozone 
problem. A simpler diurnal averaging procedure is needed so that a constant 
sun integration reproduces the true diurnal average of the long-time integration 
of a fully time-dependent calculation. In this work, the averaging procedure 
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suggested in reference 18 is adopted. This method is based on the fact that 
species having significant diurnal variations are driven mainly by photochem- 
istry. In addition, most of this variation occurs near sunrise and sunset, 
while over a major portion of the day and night the variation is much smaller 
and can be characterized by average values. Thus, the task reduces to determin- 
ing the ratio of the average nighttime to daytime concentration, 
which is found by appropriately simplifying the time-dependent photochemical 
equation (without diffusion). The parameters r^j^ are used to modify the photol- 
ysis and gas kinetic rates by taking a 24-hour average of the chemical terms in 
the conservation equation. For a bimolecular reaction between species i and j, 
a typical chemical term can be written as 

kninj 

and its 24-hour average is 

/ ^ 

1 f to /td 

to -*0 ^ \!^ ^ / 

where to = 24 hours and and represent the length of the day and 

night. However, n^"^ is related to n^*^ through the quantity r^, so that 
equation (26) reduces to 


(25) 


(26) 




XniS'lr - - 


(27) 


The quantity n^^ can be expressed in terms of the 24-hour average n^ as 
follows: 

1 fto t*^ ^ t" ^ /t<3 

^i = — J ni dt = — ni'3 + — ni" = ni<3 — + ri — 

^o 0 ^o ^ \ ^o ^o , 


= ^ini^ 

Use of this last expression in equation (27) gives 


(28) 



6i6j 


n^nj = kn^nj 


(29) 
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which has the same form as equation (25) . For a term involving photolysis, the 
24-hour average is 



Jini 


dt 




td 

- 

to 



(30) 


where it is assumed that the photolysis rate is constant over the day at its 
daytime average. The last term in equation (30) is obtained by noting that Ji, 
the 24-hour average photolysis rate, is related to its daytime average according 
ta 

to and by using equation (28) . Values of were calculated by 

^o 

a procedure analogous to that employed in reference 19. This method allows for 
the variation of solar zenith angle over the day rather than evaluating the 
photolysis rate at some mean zenith angle. 

As the derivation outlined in the previous paragraph illustrates, the key 
quantity for including diurnal effects in a constant sun model is the param- 
eter r^. For many species, this parameter can be determined immediately. 
Species with long chemical lifetimes, such as CH 4 , N 2 O, Freons, CO, HCl, HNO 3 , 
and O 3 , have little diurnal variation so that = 1 is a good approximation. 

On the other hand, species with very short chemical lifetimes, such as 0 and NO, 
decrease rapidly at sunset to essentially zero at night so that r^ « 0. For 
species with intermediate lifetimes, must be calculated by suitably sim- 

plifying the photochemical equations for nighttime chemistry. In the present 
work, this procedure is used to obtain explicit expressions for r^ for NO 2 , 
NO 3 , N 2 O 5 , C10N02r CIO, and HO 2 . All other species are assumed to have a value 
of r^ equal to either 0 or 1 depending on their reactivity. It should be 
noted that a species may exhibit some dirunal variation, which implies that r^ 
is near unity, but because it is relatively unimportant in the chemistry, the 
simpler procedure of ignoring its variation is adopted. In order to illustrate 
the type of manipulations required, the derivation of r^ for CIONO 2 is given 
in detail in the appendix. A summary of the results for all of the species, 
along with basic assumptions made in the derivations, is also given in the 
appendix i 

INPUT DATA 
Transport 

The eddy diffusion coefficient used in the present calculations is plotted 
in figure 3 along with others used recently in various models. The current 
choice of the function for this model is that recommended in reference 17. 

At the present time, it is difficult to state what the most appropriate K 2 
profile is for one-dimensional calculations. At best, K 2 must be regarded 
as a transport parameter determined, insofar as possible, by available tracer 
or vertical mixing data. 
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Photochemical Data 


The kinetic mechanism considered in the model includes 115 reactions 
involving 36 active species. The specific reactions and rate data employed are 
listed in table I, along with the data sources from which they were obtained. 
(See refs. 4, 12, 17, and 20 to 44.) Solar intensity data at the outer edge 
of the atmosphere were taken from reference 12, with the modifications to the 
ultraviolet wavelengths indicated in reference 13. The photolysis rates shown 
in table 1(b) are noontime values at 55 km for a 45^ solar zenith angle. Molec- 
ular oxygen number densities are held fixed and are taken from reference 10. 

It must be recognized that large uncertainties remain in many of the rate coef- 
ficients, and the possibility exists that important processes and species 
remain to be discovered. Thus, current sets of reactions and species are sub- 
ject to periodic review and alteration. 


Loss Rates for Condensible and Soluble Species 

To represent in a reasonable way the globally averaged vertical concen- 
tration profiles of condensible and soluble species in the troposphere, it is 
necessary to account for precipitation removal processes. Of course, the most 
important species in this category for the Earth's atmosphere is H 2 O, and the 
use of its precipitation factor in this model is presented in some detail. 
Below the tropopause, precipitation processes are primarily responsible for 
the average lifetime of a water vapor molecule in the atmosphere. As a first 
approximation in determining a value for this lifetime, a continuity equation 
for water vapor is written: 


^^H20 >^H20 

‘5h20 


(31) 


where 


%20 


is the precipitation lifetime of H 2 O and the last term represents 


a loss rate of water vapor due entirely to the condensation and subsequent 
rainout processes. Production and loss due to chemical reactions are negli- 
gible. Using the vertical eddy diffusion coefficient defined previously and 
assuming steady state, this equation becomes 


d j 

^"H2o\ 

-nK^ -1 
dz 

Vn jj 


nH20 

= 0 

‘^H20 


(32) 


where n is the total number density. Since the water vapor mole fraction may 
be defined as 
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(33) 


Xh20 = 


^H20 

n 


equation (32) may be rewritten as 


j 2 Y 

«20 d(ln n) ‘^^H20 XH2O 

+ 0 

dz2 dz dz Kz %20 


(34) 


This equation may be used to determine reasonable values of ^^20 
one-dimensional models if ^h 20' known functions of altitude - 


With 


values of n and XH2O reference 45 and values of Kg cited pre- 


viously ^ the inverse lifetime ^ deduced from a numerical solution of equa- 
tion (34) , is found to be closely represented by 


= K2(0.076z + 0.25)2 

%20 


(35) 


where Kg is in m^-s"*^ and z is in km. Since water vapor is relatively 
unreactive in the troposphere, this analysis should yield a value of 

which gives reasonable results in a coupled chemistry model. Specifically, 
for ambient conditions, the global average tropospheric water vapor profile 
from which the lifetime expression is derived should be essentially reproduced. 
When studies are performed involving water vapor injections at various alti- 
tudes, the present formulation provides a means of investigating resultant 
water vapor production, loss, and altered concentration profile in a self- 
consistent manner within the model. The inclusion of the lifetime expression 
in the model is simply accomplished by adding the term nH20 

chemical loss term Lh 20^H20* Figure 4 shows the tropospheric water vapor 

profile for the unperturbed atmosphere calculated with the model, along with 
profiles obtained from other sources. It is apparent that the present calcu- 
lation represents what may be reasonably considered a globally averaged vertical 
distribution. 

A very rough estimate of the water vapor lifetime may be obtained by 
dividing the average atmospheric water vapor content by the mean precipita- 
tion rate. The total water vapor content may be easily approximated by 


/%20 
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integration of a mole fraction profile represented by a 2 -km scale height 
distribution. (This profile is also shown in fig. 4.) Such a calculation 
leads to a value of 2.67 x 10^^ H 2 O molecules in the atmosphere. Reference 46 
gives the annual precipitation rate as 200.9 cm-year"^ , which converts to 
1 .09 X 10^8 molecules-s ^ • The average lifetime based on these two values is 


2.83 days. The tropospheric lifetime q. 


H2O 


obtained from equation (35) varies 


with altitude between values of about 1 to 1 8 days. Thus, an average lifetime 
as represented in equation (35) is reasonably similar in magnitude to the value 
based on data sources. However, it is impossible to state exactly what the 
global mean lifetime of H 2 O actually is without access to more accurate and 
widespread data on global precipitation. As is stated in reference 47, the 
atmosphere holds on the average about 10 days’ supply of water, but, in heavy 
rainfall areas, the columnar H 2 O content demands a moisture supply v^hich cor- 
responds to a lifetime of only 1 or 2 hours. 


In addition to water vapor, the soluble species for which rainout fac- 
tors are included in the troposphere are HCl, HNO 3 , and CH 3 OOH. Of course, 
for these trace constituents, an estimate of rainout lifetime as described 
previously for water vapor is not possible because of their relative chemical 
activity in the troposphere and a lack of sufficiently extensive global data 
on their average vertical concentration profiles. For these reasons, the rain- 
out factors presently used for these species are rough estimates. A constant 
value of inverse rainout lifetime of 3 x 1 0“^ s"^ is used for both HNO 3 and HCl 
(ref. 48) for 0 < z < 9 km, while the precipitation factor for CH 3 OOH is equated 
to that of H 2 O. More rigorous approaches to the problem of soluble gases in the 
atmosphere may be found, for example, in reference 49. 


Radiation Data 

Figure 5 summarizes the solar and infrared bands of O 3 , H 2 O, and CO 2 
treated in the radiative transfer model. Band absorptance and emissivity 
expressions, along with required input data such as bandwidths and strengths and 
line half -widths, are taken from reference 8 . 


Boundary Conditions 

In addition to the chemical and radiative input data described previously, 
conditions at the upper and lower boundaries of the calculations are of utmost 
importance in this type of solution; and for some of the minor atmospheric con- 
stituents, specification of appropriate boundary conditions is difficult. In 
this model, three types of boundary conditions may be used: fixed number den- 

sity, fixed particle flux, and photochemical equilibrium. The boundary con- 
ditions imposed on each of the species at 0 and 55 km for the present calcu- 
lation are given in table II. When chemical lifetimes are relatively short (on 
the order of minutes) , equilibrium conditions are used, and when number den- 
sities have been well established, these values are selected. The specification 
of zero flux for a species means that the volume concentration, or mole frac- 
tion, of that species is constant across the boundary. 


SAMPLE CALCULATIONS 


Using the input data and boundary conditions just described, some results 
of a temper ature-coupled solution of the current model are presented. Although 
these results are for a nominal ambient atmosphere, the model is capable of 
including chlorofluoromethanes (CF 2 CI 2 and CFCI 3 ) and their by-products in the 
calculations. Figure 6 shows the vertical concentration profiles of the various 
species for the converged steady-state solution. The discontinuity in the 
parameter at 14.5 km (see fig. 3) results in slope discontinuities in several of 
the species profiles. The cloud cover parameters at 5 km result in similar 
behavior. These profiles have been found in general to compare favorably with 
other computations and experimental data distributions. Figure 7 compares the 
temperature profile calculated by the model with the profile of reference 10 . 

The temperature deviations in the lower stratosphere are characteristic of one- 
dimensional radiative-convective calculations (see for example, ref. 8 ) and are 
due to dynamical effects which are not taken into account in the model. 

The program output includes printed values of number densities and/or mole 
fractions for all species at each altitude node point as well as the temperature 
structure. Optional printout and plotting is available for the production and 
loss of each species due to each reaction. Curves of this sort for the various 
species are extremely helpful in analysis of the solution, for example, in deter- 
mining which processes are significant at various altitudes. 


CONCLUDING REMARKS 

A computational one-dimensional model of the Earth's atmosphere has been 
discussed and some of the more salient and unique features have been described 
in detail. In addition, input data based on the most recent available infor- 
mation have been presented, along with results obtained using these data. This 
model, along with several others in current use, gives a reasonable representa- 
tion of many aspects of the globally averaged atmosphere and should continue to 
be useful in future studies of processes which impact the physical and chemical 
state of the atmosphere. Of course, of primary concern are those perturbations 
which impact the ozone layer in the stratosphere and the overall thermal balance 
of the planet. 

It must be recognized that a one-dimensional model, however refined, is 
only a first approximation to the actual globally averaged atmosphere. Never- 
theless, it is encouraging that this and similar models depict many important 
chemical and thermal properties in satisfactory agreement with observations. 
Because of this feature, their relative simplicity and their computational 
economy, one-dimensional models will have continued usefulness in a variety of 
applications, such as parametric and sensitivity studies of various chemical 
schemes and initial assessments of man-induced perturbations of the physical 
state of the atmosphere. 

Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 12, 1 979 
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APPENDIX 


DERIVATION OP DIURNAL PARAMETERS FOR C10N02r 
CIO, HO2, NO2, N2O5, AND NO3 

For the Langley model, the six species CIONO2, CIO, HO2, NO2, N2O5, and NO3 
are the only ones considered for special values of the diurnal parameter 
The other species are assumed to have either no diurnal variation or to be pres- 
ent only during the presence of sunlight. In the latter category are NO, 0, OH, 
H, 0(**D), N, CIO2, CH3, HCO, and CH3O. Those species which are assumed not to 
vary between day and night are O3, HNO3, HNO2, H202» N2O, H2O, HCl, CI2/ CF2Cl2f 
CFCI3, CH4, CO, CCI4, CH3CI, CH2O, and CH300H, The output of the calculations 
can be formulated to give either daytime concentrations and/or the average 
24-hour concentrations. 

In this appendix the diurnal parameter for CIONO2 is derived, and the 
derivations of r^ for the other five species are summarized. The diurnal 
parameters derived pertain specifically to the set of reactions listed in 
table I. Inclusion of additional species or changes in the rate constants 
will require a reevaluation of these parameters because the relative ordering 
of the photochemical terms retained may be changed. The procedures outlined 
in this appendix, however, remain unchanged. 


Derivation of Diurnal Parameters for CIONO2 and CIO 

The only significant nighttime chemistry for CIONO2 involves its for- 
mation due to the reaction between CIO and N02, so that its time variation 
during the night can be written approximately as 

d[C10N02l 

« k57[N02]*^[M] [CIO] = pn[ciol (A1 ) 


while during the day, its variation is approximately 


d[C10N02l 

= P<3[cio] - L«3[C10N02] 


(A2) 


where 


= k 67 [N 02 l'^[M] 


L<3 = J21 + k69[0]<3 


(A3) 
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and the superscripts d and n signify daytime and nighttime averages. In 
order to solve these equations, it is necessary to relate the concentration 
of CIO to that of C 10 N 02 * This relationship is obtained by noting that the 
total Cl family concentration, tci^l r which includes all species containing 
chlorine other than the source molecules (Freons, CCI 4 , and CH 3 CI) , is long- 
lived and has no appreciable diurnal variation. Thus, over 24 hours 


[Cl^l = [Cl] + [CIO] + [ClOO] + 2[Cl2] + [HCI] + [CIONO 2 ] (A4) 


is essentially constant. Although the CI 2 concentration increases at night, 
its contribution to the family sum is a few percent at most, so that its 
variation is neglected and its concentration is set equal to its daytime 
average value. Furthermore, exact time-dependent calculations (ref. 16) show 
that HCI has no significant diurnal variation (i.e., rnci = f)* Use of these 
last two assumptions in conjunction with equation (A4) allow one to write 

[cix]| = [Cl] + [cio] + [ClOO] + [CIONO 2 ] (A5) 

= Constant 


Equation (A5) is further simplified by assuming that [CIO] , [Cl] , and [ClOO] 
rapidly equilibrate. Thus, during the day, equation (A5) reduces to 

[cij^ = [ClOlG + [CIONO 2 I (A 6 ) 


where 


G = 1 + f 
f = fT (1 + f2) 

k45[0]«3 + k46[NO]^ 
^ k44[03] 


f2 


ks5 [O 2 ] 


•<56 


At sunset, [Cl] and [ClOO] go very rapidly to zero while [ciol changes more 
slowly, so that for a few hours after sunset. 



[cio] + [CIONO2] 


(A7) 
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Equations (Al ) and (A2) can now be solved using equations (A6) and (A7) and 
the constraint that [CIONO 2 ] at time t + 24 equals its value at time t. The 
resulting expressions can be integrated to obtain daytime and nighttime average 
concentrations. After much algebraic manipulation, the final results are 


pd 

[CIONO 2 I ^ ^ 

+ F-]L<3 


lT 

[c10N02]'^ 

L<3 

r- 1 ^ ■ 

- ^2 

[Clx] 



where 


pd 

lT = - -h l<3 
G 


[l - exp(-P’^T'^) ] [l - exp(-LTT<3)] 


lTt< 3{1 - exp[-(LTT<3 + PHtO) ] } 


lTt<3 

P2 = Fi 

pnrpn 


(A8) 


(A9) 


Daytime and nighttime average concentrations of CIO can be expressed in 
terms of those for CIONO 2 by using equations (A6) and (A7) as follows; 


[CIO] <3 _ 1 / [C10N02l<3\ 


(AlO) 


[CioJn [CIONO2]" 

[^x] [^x] 
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The ratios of equations (A9) with (A8) and (All) with (AlO) give the values 
of ri for CIONO 2 and CIO, respectively. 


Derived Diurnal Parameters for Other Species 

NO 2 .- For this species, it is assumed that the daytime value of [noI + [NO 2 ] 
is equivalent to the nighttime value of [NO 2 ] (see ref. 16); thus. 


^N02 


[n02]*^ 

[NO 2 ] ^ 


[no] <3 + [n02]^ 

= 1 

[NO 2 ] ^ 


[no] <3 

+ 

[NO 2 ] <3 


(A12) 


Using the dominant production and loss terms in a chemical equilibrium expres- 
sion for [NO 2 ] gives 


[no] <3 + kg[0]d 

[n 02]<3 ki3[03]<3 + k2otH02]<3 + k46[C10]<3 


(A13) 


which may be substituted directly into equation (AT 2) . 

HO 2 .- To obtain the diurnal parameter for this species, it is assumed that, 
during daylight, photochemical equilibrium conditions exist. During the night, 
retention of the dominant chemical terms yields the following rate expression: 


d[H02] 

— = -Ln[H02] - 2k22tH02]2 


(A14) 


where 


Ln = k2s[03] + k7o[N02]n 


Equation (A1 4) may be integrated over the nighttime t^ and an average night- 
time concentration obtained; 


1 

[H 02 ]^ = In 

2k22t^ 


2k22[H02]^ 

+ [1 

LD 


exp (-L*^t’’) ] 


[h02]<3 

[1 - exp(-L>’t*^) ] (A15) 

Lnt^ 
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This last approximation is valid because 


21^22 [HO 2 ] ^ 

« 1 

IP 


Thus, 


^H 02 


1 


[1 


Lntn 


exp ] 


(AT 6) 


N 2 O 5 .- Two differential equations for rate of change of N 2 O 5 (a daytime 
equation and a nighttime equation) are solved by using the dominant production 
and loss terms. During sunlight hours. 


d[N 205 ] 

' — S3 

dt 



kjelMl 


+ k38 [0] *3 ] [N2O5] 


(A17) 


During nighttime (assuming chemical equilibrium for [NO 3 ]) , 


d[N 205 ] 

— — IStf 

dt 


ki4[N02l"[03] 


(AT 8) 


When equation (AT 7) is integrated over daytime and equation (AT 8 ) over nighttime 
and the constraint of cyclic variation is imposed, the following expressions 
are obtained: 


ki 4 [n 02] >^[03] tn 

[N205]<3 = (A19) 


[N 2 O 5 ] n 


1 

2 


ki 4 [n02]'^[03] 


1 + exp(-L<3t<3) 
1 - exp(-L<3t^) 


(A20) 


Thus, 
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''N2O5 


1 + exp(-L‘3t‘3) 

2 1 - expC-L'^t'^) 


where 


L<a = Jt 3 + k3g[M] 


+ k38[0l ^ 


NO3 For the diurnal parameter of NO3, chemical equilibrium is 
for both conditions. In the absence of sunlight. 


[n 03]’^ « 


)<3g [N2O5 ] '^TmI + k-|4 [NO2I ’^[03] 
k35[N02l "Tm] 


The daytime expression is more complex: 

d 

%03 

[m]6 f 

d 

%03 


where 

d d 

PNO3 = Jl3rN205l<3 + kg9[0] <3[ci0N02l + kgi [cil <3 [hN03] 

+ k 3 g[M] [N 2 O 5 ] ^ + k 34 [N 02 l <3 [m] + k29 [OH] <3 [hN03] 

+ k-| 4 [ 03 l [n 02]^ + k 9 [o]^[HN 03 ] 

and 

d d d 

%03 = Jn + JlO + k35fN02l°[M] + k23[NOl<3 
As before, the ratio of equations (A22) to (A23) gives the value of 


(A21) 


invoked 


(A22) 


(A23) 


^N03* 
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TABLE I.- PHOTOCHEMICAL PROCESSES FOR LANGLEY ONE-DIMENSIONAL MODEL 


(a) Kinetic process 


Process 


( 1 ) 0 ('*D) + O2 -*■ 0 + O2 

( 2 ) 0 ('' D) + N2 -*■ O + N2 

( 3 ) 0 (^ D) + H2O ^ 20 H 

( 4 ) O + O2 + M ^ O3 + M 

( 5 ) O + O3 2O2 

(6) O + OH H + O2 

( 7 ) O + HO2 OH + O2 

(8) 0 + NO2 ^ NO + O2 

( 9 ) 0 + HNO3 -»• NO3 + OH 

( 10 ) O + HNO2 NO2 + OH 

( 11 ) O2 + H + M -»■ HO2 + M 

( 12 ) O3 + OH ->■ HO2 + O2 

(1 3 ) O3 + NO -»■ NO2 + O2 
(1 4 ) O3 + NO2 ^ NO3 + O2 
( 15 ) O3 + H OH + O2 
(1 6) 20 H -»■ H2O + O 
( 1 7 ) OH + HO2 ^ H2O + O2 

(1 8) OH + NO2 + M HNO3 + M 

( 19 ) OH + NO + M HNO2 + M 

( 20 ) HO2 + NO -»■ OH + NO2 
( 21 ) OH + H2O2 H2O + HO2 

(22) 2H02 ^ H2O2 + O2 

( 23 ) NO + NO3 2 N 02 

( 24 ) NO + 0 + M ^ NO2 + M 

( 25 ) N2O + 0 (^ 0 ) N2 + O2 

( 26 ) N2O + 0 ( 1 d) ->■ 2 NO 

( 27 ) 0 (^D) + N2 + N2O + M 

( 28 ) HO2 + O3 ^ OH + 2O2 

( 29 ) OH + HNO3 H2O + NO3 

( 30 ) 0 (^ D) + H2 OH + H 

( 31 ) OH + HNO2 H2O + NO2 

( 32 ) 0 (^D) + CH4 -»• OH + CH3 

( 33 ) NO + HO2 + M HNO3 + M 

( 34 ) NO2 + 0 + M -*■ NO3 + M 

(35) NO2 + NO3 + M ^ N2O5 + M 

( 36 ) N2O5 + M ^ NO2 + NO3 + M 

( 37 ) N2O5 + H2O 2HNO3 

( 38 ) N2O5 + O ^ 2 N 02 + O2 

( 39 ) N + NO2 N2O + O 

( 40 ) N + O2 -»■ NO + O 

See footnotes at end of table. 


Rate expression 

Ref. 

(a) 


2.9 X lO-ll exp(67/T) 

20 

2.0 X 1 O"''"' exp(l 07/T) 

20 

2.3 X 10-‘*0 

20 

1.07 X 10~34 exp(510/T) 

20 

1.9 X 10“''"’ exp(-2300/T) 

20 

4.2 X 10-’’^ 

20 

3.5 X lO-ll 

20 

9.1 X 10-'' 2 

20 

2.0 X 1 0-'>7 

20 

U) 

• 

0 

X 

0 

20 

2.08 X 10'’32 exp(290/T) 

20 

1.5 X 1 0-“' 2 exp(-l 000/T) 

20 

2.1 X 1 0 ”^ 2 exp{-1 450/T) 

20 

1.2 X 10-’’3 exp(-2450/T) 

20 

1.0 X 10“^^ exp(-516/T) 

20 

1.0 X 10”^^ exp(-550/T) 

20 

3.0 X lO-ll 

20 

t>0.94 <1)1 ([m1 ,T) 

20 

0.25kiB 

(c) 

8.1 X 10-' 2 

20 

1.0 X 10“^^ exp(-750/T) 

20 

2.5 X 1 0""’ 2 

20 

1.9 X 1 O-ll 

20 

1.55 X 1 0“32 exp(584/T) 

20 

5.5 X 1 O-l^ 

20 

5.5 X 1 O-l'' 

20 

3.5 X 1 0-27 

20 

7.3 X lO-l^ exp(-1 275/T) 

20 

r- 

1 

0 

X 

0 

00 

20 

9.9 X 10“*^' 

20 

6.6 X 10-’’ 2 

20 

1.3 X 1 O-l 0 

20 

1.0 X 10-'' 5/[m1 

20 

1.0 X 10-3'' 

20 

^< 1)3 ([m1 ,T) 

20 

®<1>4([m] ,T) 

20 

1.0x1 0-20 

20 

1.0 X 1 0-'' 4 

20 

2.0 X 10-"''' exp(-800/T) 

20 

5.5 ^ 10-"'2 exp(-3200/T) 

20 
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TABLE I.- Continued 


(a) Continued 


Process 

Rate expression 
(a) 

Ref. 

(41) N + NO ^ N 2 + 0 

8.2 X 10-'*'' exp(-410/T) 

20 

(42) N + NO 2 ^ 2N0 

0 

20 

(43) N + O 3 ^ NO + O 2 

2.0 ^ lO-l^ exp(-1070/T) 

20 

(44) Cl + O 3 ^ CIO + O 2 

2.7 X 10"'''' exp(-257/T) 

20 

(45) CIO + 0 -► Cl + O 2 

7.7 X 10"'' 1 exp(-130/T) 

20 

(46) CIO + NO -»• Cl + NO 2 

1.0 X 10"^^ exp(200/T) 

20 

(47) Cl + CH 4 -»■ HCl + CH 3 

7.3 X 10 "'' 2 exp(-1269/T) 

20 

(48) Cl + H 2 HCl + H 

3.5 X 10"'''' exp(-2290/T) 

20 

(49) Cl + HO 2 -> HCl + O 2 

3.0 X 10 "'''' 

20 

(50) Cl + H 2 O 2 -»■ HCl + HO 2 

1.7 X 10 "^ 2 exp(-384/T) 

20 

(51) Cl + HNO 3 ■> HCl + NO 3 

1.0 X 10"^^ exp(-2170/T) 

20 

(52) H + HCl ^ H 2 + Cl 

7.8 X 10"'' 2 exp(-1600/T) 

20 

(53) HCl + (M ^ Cl + H 2 O 

3.0 X 10"'' 2 exp(-425/T) 

20 

( 54) HCl + 0 ^ Cl + OH 

1.14 X lO"''! exp(-3370/T) 

20 

(55) Cl + O 2 + M •> CIO 2 + M 

1.7x1 0"33 

20 

(56) CIO 2 + M Cl + O 2 M 

5.8 X 10"09 exp(-3580/T) 

21 

(57) 2C1 + M ^ CI 2 + M 

6.0 X 10"34 exp(900/T) 

20 

(58) Cl + CIO 2 2C10 

1 .1 X 10 "'''' 

20 

(59) Cl + CIO 2 ^ CI 2 + O 2 

• 

C\ 

X 

0 

1 

c 

20 

(60) 2C10 + M ^ CI 2 + O 2 + M 

0 

17 

(61) CI 2 + 0 CIO + Cl 

4.2 X 10"'' 2 exp(-1370/T) 

20 

(62) Cl + CIO CI 2 + 0 

0 

(f) 

(63) 2C10 > CIO 2 + Cl 

2.0 X 10 "^ 2 exp(-2200/T) 

20 

(64) 2C« + M ^ H 2 O 2 + M 

1.25 X 10"32 exp(900/T) 

20 

(65) H 2 O 2 + 0 •> CH + HO 2 

2.75 X 10 "'' 2 exp(-2125/T) 

20 

( 66 ) C«J + CH 4 -»■ CH 3 + H 2 O 

2.36 X 10"'' 2 exp(-1710/T) 

20 

(67) CIO + NO 2 + M ^ CIONO 2 + M 

9(1)5 ([M] ,T) 

20 

( 68 ) HCl + CIONO 2 CI 2 + HNO 3 

0 

3.0 X 10 "^ 2 exp(-808/T) 

(f) 

(69) C 1 C »102 + 0 CIO + NO 3 

20 

(70) HOo + NO 2 HNO 2 + O 2 

(71) 0('D) + HCl ^ OH + Cl 

2.0 X 10"15 

20 

• 

X 

0 

1 

c 

20 

(72) H 2 + OH -»■ H 2 O + H 

3.6 X 10"** exp(-2590/T) 

20 

(73) 2C10 CI 2 + 02 

5.0 X 10"'' 3 exp(-1238/T) 

20 

^(74) 0 ("'d) + CF 2 CI 2 ^ 2C1 + . . . 

2.0 X 10"'' 0 

20 

h(75) 0(^D) + CFCI 3 2C1 + . . . 

2.3 X 10"'' 0 

20 

(76) CH 3 + 02-^ CH 2 O + OH 

2.9 X 10 "! 3 exp(-940/T) 

20 

(77) CH 3 + O 2 + M ^ CH 3 O 2 + M 

1.9x1 0"31 

20 

(78) CH 3 O 2 + HO 2 CH 3 OOH + O 2 

6.7 X 

20 

(79) CH 3 O 2 + NO -»■ CH 3 O + NO 2 

3.3 X 10"! 2 exp(-500/T) 

20 

(80) CH 3 O 2 + NO -»• CH 2 O + HNO 2 

0 

20 


See footnotes at end of table, page 33, 
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TABLE I.- Continued 


(a) Concluded 


Process 

Rate expression 
(a) 

Ref. 

(81) CH 3 O 2 + NO 2 CH 2 O + HNO 3 

0.25 (R79/2.2) 

22 

(82) CH 3 O + 02 ^ CH 2 O + HO 2 

1.6 X 10-''2 e3^(-3300/T) 

20 

(83) CH 2 O + OH -»• HCO + H 2 O 

3.0 X lO-l^ exp(-250/T) 

20 

(84) HOO + O 2 ^ 00 + HO 2 

5.7 X 10-12 

20 

(85) CO + OH -> CO 2 + H 

i(j, 6 (fM]) 

23 

( 86 ) CH 3 CI + OH ->• CH 2 + Cl + H 2 O 

2.2 X 10-12 exp(-1142/T) 

20 

(87) CH 3 OOH + OH -»■ CH 3 O 2 + H 2 O 

k 21 

24 


®Bimolecular rates are in units cmVparticle-s, termolecular rates 
in cm®/particle^-s. 

^log-|o 4>i = -AT/(B + T) - 0.5 logiQ(T/280) - ? where 
C = log 1 0 [m] 

A = 31.62273 - 0.258304? - 0.0889287?2 + 0 . 0025201 73?^ 

B = -327.372 + 44.5586? - 1 .38092?2 

‘^This expressicxi closely aEproximates tabular rate functicxi given in 
reference 20. 

= (1/[m1(1.48 X 10 “3) exp(861/T) 

® 4>4 = 4>3/(1.27 X 10-27 exp(11180/T) 

^Rates which are unimportant and/or for which insufficient 
information is available are equated to zero. 

9 <j >5 = (3.3 X 10-23 T-3-34)/(t + s.7 x 10-^ T-0*6 [m]V2) 

^The only active products considered for these processes are 
chlorine atoms. 

i<t >6 = (1.5 X 10~‘'2 + 9.316 x 10-33 [m])/( 1 + 1 .5526 x 10-20 [m1) 
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TABLE I.- Continued 


(b) Photolytic processes 


Process 


( 1 ) O 2 + hv 20 

(2) O3 + hV -► 0 + O2 

O3 + hv 0(^D) + O2 
NO2 + h\) NO + 0 
NO2 + hV -> NO + 0(^D) 
HNO3 + hV -► NO2 + OH 
HNO2 + hV NO + OH 
H2O + hV -► H + OH 
(9) H2O2 + hV 20H 

(10) NO3 + hV NO + O2 

(11) NO3 + hV ->• NOo + 0 

(1 2) N2O + hV 0(’D) + N2 


(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 


03) 

(14) 

•05) 


( 21 ) 

( 22 ) 


Wavelength 
range, nm 


Computed values of J, s”^ , at z, km, of 


N2O5 + hv NO2 + NO3 
NO + hv -► N + 0 
HCl + hV H + Cl 

(16) CIO2 + hv •> Cl + O2 

(17) CIO + hV Cl + 0 

(18) CI2 + hv 2C1 

(19) CFCI3 + hv -> 2.5C1 + . . 

(20) CF2CI2 + hV -»* 2C1 + . . 

CIONO2 + hv -► CIO + NO2 

CCI4 + hv 2C1 + . . . 

(23) CH3CI + hv -► CH3 

(24) CH2O + hv •+ 

(25) CH2O + hV 


> CH3 + Cl 
H + HCO 
CO + H2 


(26) 

(27) 


HCO + hV CO + H 


CH3O2 


+ hv CH^O + 0 


(28) CH3OOH + hv CH3O + 


OH 


3116. 3 
116.3 
116.3 
130.7 

130.7 
190.5 

303.0 

116.3 

175.4 

397.5 

397.5 

116.3 

210.5 

183.0 

138.9 

224.7 

227.3 

246.9 
185.2 
185.2 
185.2 

I 185.2 
' 173.9 

289.9 
289.9 

175.4 
175.4 


to 235.3 0. 
to 735.0 
to 735.0 
to 497.5 
to 497.5 
to 322.5 
to 397.5 
to 198.0 
to 347.5 
to 705.0 
to 705.0 
to 250.0 
to 377.5 
to 191 .0 
to 219.8 
to 281 .7 
to 281.7 
to 447.5 
to 227.3 
to 222.2 
to 462.5 
to 227.3 
to 219.8 
to 362.5 
to 362.5 
(b) 

to 347.5 
to 347.5 


294 X 

380 X 

>154 X 

,110 X 

,641 X 

.392 X 

.607 > 

.654 > 

.783 > 

.261 > 

.261 > 

.939 > 

.292 > 

,549 

.615 

.491 

.929 

.246 

.136 

.168 

.637 

.173 

.262 

.320 

.900 

.261 

.783 


10-23 
10“3 
10-5 
10-1 
10-19 
10“6 
10-5 
10-22 
10-4 
10-2 
10-1 
10-20 
10-4 
c 10-48 
c 10-20 
( 10-28 
( 10-28 
< 10-2 
< 10-1® 

< 10-19 

< 10-4 

< 10-18 
< 10-20 
X 10-4 

X 10-4 

X 10-4 
X 10-4 


0.122 

.388 

.195 


10-19 

10-3 


16 


0.587 

.398 

.890 


10-15 

10"5 

10-5 


25 


0.233 

.428 

382 


10-11 

10-5 

10-4 


40 


0.519 

.461 

.183 


10-9 

10-3 

10-2 


55 


Refs. 


o.no 

.460 

.822 


10-8 12, 24, 25 
10-3 12, 26 

10-2 12, 26 


.110 

X 

10-^ 

.110 

X 

10"’' 

.112 

X 

10"'' ' 

.116 

X 

10"'' . 

.117 

X 

10"'' 

27 

.256 

X 

10-15 

. 1 21 

X 

10-10 

.483 

X 

10-'7 

.157 

X 

10-4 

.347 

X 

10"4 

28 

.412 

X 

10"® 

.467 

X 

10"® 

.142 

X 

10-5 

.569 

X 

10-4 

.990 

X 

10’4j 

29 

.607 

X 

10-3 

. 609 

X 

10“3 • 

,618 

X 

10“5 

.632 

X 

10"5 

.633 

X 

10-^ 

29 

.404 

X 

10-19 

.466 

X 

10-15 

.442 

X 

10-11 

.209 

X 

10-8 

.141 

X 

10"7 

27 

.791 

X 

10-4 

.813 

X 

10-4 

.944 

X 

10-4 

.213 

X 

10-5 

,372 

X 

10"5 

30, 31 

.262 

X 

10-2 

.268 

X 

10-2 

.273 

X 

1 0-2 

.273 

X 

10-2 

.273 

X 

10"2 

29 

.262 

X 

10-1 

.263 

X 

10"'' 

.268 

X 

10-'' 

.273 

X 

10-'' 

.273 

X 

10"'' 

29 

.1 1 1 

X 

10-15 

.1 12 

X 

10-11 

.408 

X 

10-8 

.375 

X 

10-8 

.632 

X 

10"8 

32, 33, ; 

.297 

X 

10-4 

.306 

X 

10-4 

.371 

X 

10"4 

.210 

X 

10"5 

.495 

X 

10"5 

29 

.550 

X 

10-37 

.157 

X 

10-26 

.460 

X 

10-13 

.406 

X 

10-6 

.395 

X 

10"5 

35 

.661 

X 

IQ-l 6 

.276 

X 

10-12 

.251 

X 

10-8 

.251 

X 

10-8 

.469 

X 

10"8 

36 

.107 

X 

10-26 

.128 

X 

10-22 

.183 

X 

10-1 5 

.558 

X 

10-5 

.567 

X 

10"2 

37 

.139 

X 

10-26 

.252 

X 

10-24 

,224 

X 

10-14 

.474 

X 

10-5 

.381 

X 

10"2 

37 

.248 

X 

10-2 

.250 

X 

10-2 

.262 

X 

10-2 

.291 

X 

10-2 

.296 

X 

10"2 

37, 38 

.165 

X 

10-14 

.178 

X 

10-10 

.642 

X 

10-7 

.605 

X 

10"5 

.993 

X 

10"5 

39 

.156 

X 

10-15 

.147 

X 

10-11 

.564 

X 

10-8 

.621 

X 

10"8 

.120 

X 

10"5' 

39 

.644 

X 

10-4 

.658 

X 

10-4 

.762 

X 

10-4 

.395 

X 

10-5 

.931 

X 

10"3 

40 

,224 

X 

10-14 

.294 

X 

10-10 

.111 

X 

10-6 

.158 

X 

10-4 

.256 

X 

10"4 

41 

.257 

X 

10-16 

.250 

X 

10-12 

.944 

X 

10-9 

.968 

X 

10-7 

.198 

X 

10-8 

42 

.325 

X 

10-4 

.337 

X 

10-4 

.412 

X 

10-4 

.730 

X 

10"4 

.784 

X 

10"4 

43 

.908 

0 

X 

10-4 

.922 

0 

X 

10-4 

.101 

0 

X 

10-5 

.122 

0 

X 

10“5 

.125 

0 

X 

10"5 

43 

.264 

X 

10-4 

.271 

X 

10-4 

.315 

X 

10-4 

.710 

X 

10-4 

.240 

X 

10-5 

C4 

' .792 

1 

X 

10-4 

.813 

J 

X 

10-4 

.944 

X 

10-4 

.213 

/ 

X 

10-5 

.372 

X 

10"5 

d44 


34 


^ross sections of reference 12 are used for 116.3 = X = 175.4 nm; correlations of reference 24 are 
cross sections of reference 25 are used for 204.1 < X ^ 235.3 nm. 

^This process omitted from present calculations. 

CJ27 = J9/3 as assumed in reference 4. 

<3j28 = J 9 as assumed in reference 44. 


used for 175.4 < X ^ 204.1 nm, and 



TABLE II,- SPECIES AND BOUNDARY CONDITIONS FOR ONE-DIMENSIONAL MODEL 


[p.E. denotes photochemical equilibr iumJ 


=* 


■ 


n 


Lower boundary 

Upper boundary 

Species 






Number density. 

Flux, , 

Number density. 

Flux, 0 , 


n 

, particles-cm“^ 

particles-m 2 “s^^ 

n, particles-m '3 

particles-m^-s"^ 


6.35 X 



0 

0 ( 3 p) 

P.E. 


P.E. 


NO 

P.E. 



0 

N02 

1 .00 X 10^ 


P.E. 


HNO3 


0 


0 

HN02 


0 

p.E. 


NO3 

P.E. 


P.E. 


H2O2 


0 


0 

OH 

P.E. 


P.E. 


HO2 


0 

P.E. 


N20 

8.02 X I 0 l 2 



0 

N2O5 

P.E. 



0 

H2O 

3.34 X IqI '7 



0 

Cl 

P.E. 


P.E. 


CIO 

P.E. 


P.E. 


HCl 

2.67 X 10 ^° 



0 

CF2C12 

^0 


ao 


CFC13 

ao 


^0 


CH4 

3.48 X I 0 l 3 



0 

C1N03 


0 

P.E. 


CH -0 


0 


0 

co2 

CM 

0 

X 

0 

IT) 

CM 



0 

CH3OOH 


0 


0 

CCI4 

2.55 X io 3 



0 

CH3CI 

1 .91 X 10^0 



0 

H 






oOd) 






N 






CI2 


\ These species are computed assuming photochemical 


CH3 


equilibrium at all altitudes 



HCO 






CH3O2 






CH3 






C102 






02 


< These species have fixed vertical 

profiles 


H2 







^For the ambient atmosphere, the chlorof luoromethanes are not normally included. 
For injection of these species, flux boundary conditions are imposed. 
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Water vapor mole fraction 


ll■l■lll■l|LJ!l■llllll■UJ III 


III II I II I II II 


□ From ref. 11 

o From ref. 45 


2 km scale height profile (see text) 

Present model calculation 



Figure 4.- Vertical profiles of water vapor mole fraction. 
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Figure 5.- Infrared and solar absorption bands used in one-dimensional 

transfer calculation. 
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Figure 6.- Vertical profiles of minor atmospheric constituents as computed 
by the Langley photochemical model, n^ is in particles-cm“^. 




rnrmiTnfrii 



Altitude, km 






Altitude, km 


X 10‘ 



logloCn^) 


OJ 


(c) Chlorine-containing species. 
Figure 6.- Continued. 
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Figure 7.- Vertical temperature distribution (dashed curve) obtained 
from coupled solution compared with U.S. Standard Atmosphere/ 

1962 (ref. 10) profile (solid curve). 
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